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TECHNICAL NOTE NO . 1193 

EFFECT OF EEODUCT OF INERTIA ON LATERAL STABILITY 
By Leonard Sternf leld 


SUMMARY 


A theoretical investigation was made to determine the effect 
of the product-of-inertia terms in the lateral equations of motion 
on the lateral-stability boundaries* The product of inertia 
results from the inclination of the principal longitudinal axis of 
the airplane to the flight path* 

The results of the calculations indicated that the product- 
of-inertia terms should be included in the lateral equations of 
motioxi to determine the lateral stability of an airplane* The 
value of the directional- stability derivative C„ required for 

(3 

stability, as determined from the calculations which Include the 
product of inertia caused by the inclination of the principal 
longitudinal axis above the flight path, is considerably less than 
the value predicted by calculations neglecting the product-of- 
inertia teamisa ujkdn p>cJ r i . 

INTROLUCTICN 


Recent fll^t tests and tests in the Langley free-fli^t tunnel 
indicated a discrepancy between the observed and calculated conditions 
for dynamic lateral stability. It was suggested that -tlie discrepancy 
mi^t be the result of neglecting terms jja the equations of motion 
invo].vlng the product of inertia* 

The theoretical investigations of the lateral stability of air- 
planes undertaken in the past have for the most part neglected the 
effect of the inclination of the principal longitudinal axis of the 
airplane to the flight path on the oscillatory-stability boundary* 

(See references 1 to 3») In the appendix of reference 1 , Zimmersnan 
mentions that the results of supplementary calculations indicated 
that to neglect the angularity of the principal axis to the fll^t 
path did not seriously affect the oscillatory- stability boundary in 
the normal-flight range and gave slightly conservative resxilts. The 
results of reference 1, consequently, are based on lateral equations 
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of motion from which the effect of the incllnatlcm of the principal 
azls to the flight path are omitted, namely, the terms Including the 
product of inertia. 

For the range of airplane parameters Investigated hy Zimmerman, 
the product- of “inertia terms probably have a small stabilizing 
effect on the oscillatory- stability botindary. However, the effect 
of product of Inertia on the results of stability calcilations 
covering the range of parameters of present-day airplanes designed 
for high-speed, hi^-altitude flight - that is, hl^ values of the 
relative-density factor u, a large increase in the effectlve- 
diheiral parameter C, caused by the use of swept-back wings, and 

the change in the mass distribution of the airplane - has not been 
investigated. 

In the present investigation, calculations were made to determine 
the effect of the product- of inertia on the requirements for lateral 
stability. The results of the computations were plotted as a function 
of the directional-stability derivative C„ and the effective- 

dihedral parameter C. . 


SYMBOLS 


V 

P 

q 

b 

s 

w 

m 





airspeed, feet per second 

mass density of air, sltigs per cubic foot 

dynamic pressure, poxmds per square foot 

wing span, feet 
wing area, square feet 
weight of airplane, pounds 
mass, slugs (w/g) 

acceleration of gravity, feet per second per second 

relative-density factor / ^ N 

\pSbJ 

radius of gyration about principal longitudinal axis, feet 
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k„ radiTis of gyration alDout principal normal axis, feet 
^o • • . 

moment-of-inertia coefficient aliout principal longitudinal 
■^o 

axis 

moment- of -Inertia coefficient about principal normal axis 
^o 

1-^ moment-of-inertia coefficient about fli^t-path axis 


Ig- moment-of-inertia coefficient about axis normal to flight 
path 

%Z pro6.uct-of- inertia coefficient vith respect to flight path 
and axis normal to flight path 


Cl 


P 

r, t 


lift coefficient 


/ JL cos y 

\qs 


rolling-moment coefficient 


/ Polling momentN 

\ qSb ) 


yawing-moment coefficient /Tawing moment^ 

V qSb 

lateral-force coefficient 

' \ gs 


angle of bank, radians 
azimuth angl.e, radians 
angle of sideslip, radians 

yawing angular velocity, radians per second (dt/dt) 
rolling angular velocity, radians per second (d^/dt) 



effective-dihedral derivative, rate of change of 

rolling-moment coefficient with angle of sideslip. 


per radian 


directional-stability derivative, rate of change of 
yawing-moment coefficient with angle of sideslip, 
per radian 
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C 

V 




t 

D 

A 


lateral-force derivative, rate of change of lateral-force v 
coefficient with angle of sideslip, per radian t9Cy/9w 

damping- in-yaw derivative, rate of change of yawing-moment 
coefficient with yawlng-angular-veloclty factor, 
per radian 




damping- in-roll derivative, rate of change of rolling-moment 
coefficient witn rolling-angnlar-velocity factor, 


per radian 




rate of change of yawing-mcaaent coefficient with rolling- 
angular -velocity factor, per radian 

rate of change of rolling-moment coefficient with yawing- 
angular- velocity factor, per radian 

rate of cMnge of lateral -force coefficient with rolling- 
angular- velocity factor, per radian 


rate of change of lateral -force coefficient with yawing- 
angular-velocity factor, per radian 

time, seconda 

differential operator (d/dt) 
angle of sweephack, degrees 


Tj angle of attack of principal longitudinal axis of airplane, 

positive when principal axis is ahove flight path, degrees 
(see fig. 1) 

7 angle of fli^t path to horizontal, positive in a climh, 
degrees (see fig. l) 

6 angle "between reference axis and horizontal, positive when 
i’eference axis is above horlzontail, degrees (see fig. l) 
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6 angle l)etween reference axis and principal axis^ positive 

when reference axis is ahove principal axis, degrees 
(see fig. 1) 

R Routh ' s discriminant 


EQUATIONS OF MOTION 


The linearized equations of motion, referred to the axes in 
figure 1, used to calculate the spiral- and oscillatory-stahility 
■boundaries for any flight condition, are! 


Rolling 


Ix 


'"l^ COS^Tl + I„ 
. -^o ^ 


sin^il) 
o / 


Tf - 


C, D 






' ^2 


- I^ ) sin T} cos T) U + D 


t - C p = 0 

H 


Yawing 




"Irj cos^'n + I^ sin^^l^D^ "0, . D 


n 




t 






/l - Ijj. jsin n cos n D 

Vo 0/ ■ ij) 


- C_ p « 0 


Sideslipping 


i Dp) - + Cy^ D^9^ -(^ tan 7 + Cy^ D^ t “ Cy^P = 0 
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■wheqre 


and 




Tho aymbols and I™ reprssonfc the momants of Inertia atout the 

Aq tlQ 

principal axes of the airplane. If the longitudinal principal axis is 
inclined at an angle ri to the flight path, the moments of inertia 
ahovit the flight-path axis and the axis normal to the flight path are: 

= Ij-^ coe^Tj + I2 sin^ 


" ^Zq nos^ + sin^ 
and the product of inertia is ’I ' 


/x^ - X, J 


/ rr f 

(. 


: _J.rri-r,)+0" i'7 
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When .9^06^^ is substituted, for it^Qe^'^ for and 3©®^^ 

for 3 in the equations written in determinan.t fom, X must he a 
root of the equation 


AX,^*- + + CX^ + EX + F = 0 


where 






- I^, §M ■ ^X^Z^Yp 


C = ^XZ^YpC? j + Ixz ^ + Jxz%^xi^ + ^X^TpCn^ 


+ Ix ^ Cnp + ^ ^ 

® “ "%Z%^L ■ Iz^ZpCL “ " ^20% ^ 
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The conditions necessary to obtain the oscillatory-stability 
boundary are that the coefficients A, B, C, and E must be positive 
and Routh's discriminant, R s= BCE - A# - B%, must equal zerot 
The spiral-stability boundary is foimd by setting P = 0. The 
completely stable region is therefore bounded by the ciirves R = 0 
and F 0, which are plotted as a function of the directional- 
stability derivative C„ and the effective -dihedral derivative C, • 


Stability Derivatives and Mass Characteristics 

Calculations were made showing tiie effect of the product-of- 
inertia terms on the oscillatory-stability boundary for a hypothetical 
supersonic fighter airplane, an experimental fighter airplane, and 
a model flora in the Langley free-fli^t tunnel. The values of the 
stability derivatives and mass characteristics of the two airplanes 
and the model are given in table 1= The contribution of the tall 
to the derivatives C_ and Cj (see reference 2) was included 

P ^ 

in the characteristics of the experimental filter airplane and the 
free-fli^t- tunnel model but was neglected in the calculations of 
the hypo'ttietical supersonic filter airplane. 


RESULTS MD DISCUSSICB 


The results of the Investigation are presented in a series of 
figures which show the oscillatory- and spiral-stability boundaries 
as a function of C^^ and C, . The solid E » 0 curve of figure 2 




‘3 


represents the oscillatory-stability boundary for landing flight . 
for a hjiiothetical supersonic fightor airplane which has its 
principal axis or fuselage inclined 5° above the flight path but 
the product of inertia is assumed to be zero. The dashed 


curve in this figure is the R = 0 boundary for the same condition 
with the product- of- inertia terms included in the equations of 
motion. A comparison of both curves shows a large stabilizing shift 
in the oscillatory-stability boundary for the case in which the product 
of inertia is taken into account. As C, is increased, the value 

of required for oscillatory stability as determined by the 

calculations including l yy is considerably less than the value 
determined by the calculations neglecting The spiral- 

stability boundary plotted in figure 2 and in subsequent figures 
applies to both sets of ' calculations since this boundary is not a 



KACA TE Uo . 1193 


9 


fmctlon of the product of inertia. This result does not mean, 
hoTOver, that the rate of divergence or convergence of the spiral 
motion will he the same inasmuch as the stability equation is 
different in each case . The effect of on the E •■= 0 

boundary for the saiae airplane in cruising flight is presented in 
figure 3» The solid curve represents the case of the principal 
longitudinal axis inclined either above or below the flight path 
but neglects the effect of Ijg. The other two E = 0 boundaries 
that were calculated included the effect of prodixct of inertia. 

For the case in which the principal axis is inclined above the 
flight path, « -0.00178, the stable region is Increased; 

whereas for the case in which the principal axis is inclined below 
tlie flight path, - O.OOI78, the stable region is reduced. 


The experimental fighter airplane and the free-flight-tumel 
model were tested in flight with the principal axis inclined l6° 
and 10°, respectively, above the flight path. The 3. Itlal calcu- 
lations made on the assumption that = 0 (the solid curves in 

figs . k and 5) indicated that the airplane and model would be 
unstable in flight for the combination of 0., and C^ denoted by 

the circled points in figures 4 and 5* The flight-test results, 
however, showed that the experimental fighter airplane was stable 
and the free -flight- tunnel model was marginally stable. Subsequent 
calculations which included the product- of -inertia terms indicated 
that the E = 0 boundary for the experimental fighter airplane 
increased the stable region to include the combination of C_ 


and C, 


tested in flight and that -the E = 0 boundary for the 


free-fllght-tunnel model was shifted very close to the marginally 
stable test point. 


In general, the inclination of the principal longitudinal axis 
above the flight path causes a stabilizing shift in the oscillatory- 
stability boundary but the extent of the shift is a function of 
other airplane parameters which are still to be investigated. 


The’ solid curves of figures 2 to 5 he considered to 
represent the E = 0 boundary for the actual flight conditions of 
an airplane which has its principal longitudinal axis in line with 
the flight path (ixz = O) provided the wings are set at an angle 

of Incidence to the fuselage to obtain the lift coefficient desired 
for flight. The dashed curve on each figure would, therefore, 
represent the E = 0 boundary for the same flight conditions as the 
solid curve but with the wings set at a different angle of incidence 
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since the principal axis is inclined to the flight path. If the 
principal axis is Inclined above the fli^t path, the angle of 
incidence which corresponds to the dashed curve is smaller than the 
angle of incidence corresponding to the solid curve. If the 
principal axis is inclined below the fli^t path, however, the wing 
incidence is ].arger "than the wing incidence of the airplane in which 
the principal axis coincides with the fli^t path. For example, 
the solid curve of figure 2 is the E = 0 boundary for an airplane 
with the wings set at an angle of 15° to the fuselage and the 
principal axis in line with the flight pathj whereas the dashed 
carve is the E = 0 boundary for the same airplane with the wings 
set at a 10° angle of incidence and the principal axis inclined 5° 
above the flight path. In figure 3^ the wing incidence is 7° for 
the solid curve, 5° for the case of the principal axis above the 
flight path, and 9° for the case of the principal axis below the 
fll^t path. The solid curves of figures it and 5 are the E = 0 
boundaries for the airplane and model designed with a wing Incidence 
of l6° and 10°, respectively; for both cases the dasl 'd curve 
represents an alx-plane with a wing incidence of 0°. This interpre- 
tation applied to the solid curves of figures 4 and 5 is only 
approximately true since the angle of attack at the tail is zero if 
the wing is at l6° or 10° Incidence and the principal axis coincides 
with the line of fli^t; whereas the calculations include the effect 
of the tall at an angle of attack of l6° and 10° on the stability 
derivatives C and C, . 

^ ‘'r 

A comparison of the solid and dashed curves in each figure 
clearly indicates the increase in the oscillatory- stability region 
for an. airplane designed with a wing set at 0° incidence, thereby 
necessitating the inclination of the principal axis above the 
flight path to obtain the desired lift coefficient. 


CONCLUDING EEMAEKB 

The results of the analysis made to investigate the effect of the 
product of inertia on the lateral-stability boundaries emphasize the 
necessity of including the product -of— inertia terms in the lateral 
equations of motion to determine the lateral stability of an airplane. 
The calculations Indicate that the inclination of the principal 
longitudinal axis above the fli^t path causes a stabilizing shift in 
the oscillatory— stability boundary. 

Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ta., August 13 , 1946 
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TABEE I 

STABILirr KESIVATIVES MD MASS CHARACTEEISTICS USED 3 H STAECLITT CALCDLATIOBS 



Hypothetical auperaoiilc 
fighter airplane 

E:^rinsntal fighter airplani 

Free-fllght-tunnel model 

CondltloQ 

Landing 

Orulalng 

landing 

Landing 

W/S, ID/sq ft 

80 

80 

35 

1.85 

b, ft 

20 

20 

34 

33 

0, slugs /cu ft 

0.00238 

0.0002 

0.00238 

0.00238 

V, ft/sse 

2 €h 

1465 

173 

44 

Cl 

1 

0.372 

1 

0.8 

Pi 

5 >» 

620 

13 

6.33 

^Xol ft 

2.02 

2.02 

5.28 

0.51 

V 

9.64 

9.64 

8.03 

1.29 

Oj^, per radian 

-0.197 

-0.197 

-0.3 

-0.17 

, per radian 
r 

0.25 

0.0929 

0.25 - 0 . 04 c_ , 

P(tall) 

0.142 - 0 . 099 C„ 

P(tall) 

Cjj^ , per radian 
P 

- 0.0198 

-0.00732 

-O.CS 9 - 0 . 04 C„ 

“P(tail) 

- 0 . 04 o 6 - 0.099C_ 

P(tall) 

radian 

r 

- 1 . 47 C_ , 

”p(tall) 

- 1 . 47 C„ 

“P(tail) 

-0.086 - 0.858a, 

^P(tall) 

-0.0131 - 1.2C 

“pCtall) 

Cy , per radian 
P 

0 

0 

0 

0 

Cy t psr radian 
^r 

0 

0 

0 

0 

Cy , per radian 

-l- 33 Cn , , 

P(tall) 

■^•^"^".(tall) 

.0.4, - 4.34=^, „„ 

-0.0074 - 1.76c„ 

°P(tail) 

^lJ(fuselaga), 
per radian 

-0.25 

-0.25 

- 0.02 

0 

r, aeg 

0 

0 

0 

-9 

A, deg 

60 

60 

35 

42 
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Fig. 1 
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\ 

X 

Principal axi^ 
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Fipure /. — Sysfem of axes one/ angular 
relationships m flight. Arrows indicate 
positive direction of angles, Q-T-f 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


D ) re c ftona A ^ fabi/i /y den vaft ve jC^ 



Effective - chhea/rar/ c/erivofive 



NATIONI^L ADVISORY 
COHMITTEE FOI AERONAUTICS 


Figure Lcfferal-sfabilify boundaries for landing flight for 
hypothetical supersonic fighter airplane. 
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Direaf/onai-sfctb/Jify a/envi 



Figure 3. ~ Laieral-sfahilify boundaries for cruising fhghf for 
hypothetical superson/c f/phfer atrp/ane . 



Di re c fion aj- sfahi h fy den vafi ve j 



0 'M 708 r/£ 7f6 tZO tZ4 7Z8 -.32 t30 

Fffecfi I'e - dihedral den va ft i/e j 

F/^ure 4.— Lateral -sfcrjbilify boundaries for landing flight 
for expert mental fiyhter airplane . 
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Fjgare JT. — Lateral- 'Sfabihiy hounc/anes for landing 
flight for 'free- flight -tun net model. 


NACA TN No. 1193 Fig. 5 



